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ABSTRACT 

Spectral energy distributions (SEDs) and radio continuum spectra are presented for a northern sample of 104 extragalactic radio sources, based 
on the Planck Early Release Compact Source Catalogue (ERCSC) and simultaneous multifrequency data. The nine Planck frequencies, from 30 
\ I 1 to 857 GHz, are complemented by a set of simultaneous observations ranging from radio to gamma-rays. This is the first extensive frequency 
coverage in the radio and millimetre domains for an essentially complete sample of extragalactic radio sources, and it shows how the individual 
t-h . shocks, each in their own phase of development, shape the radio spectra as they move in the relativistic jet. The SEDs presented in this paper were 
^ ■ fitted with second and third degree polynomials to estimate the frequencies of the synchrotron and inverse Compton (IC) peaks, and the spectral 
f-H indices of low and high frequency radio data, including the Planck ERCSC data, were calculated. SED modelling methods are discussed, with an 
emphasis on proper, physical modelling of the synchrotron bump using multiple components. Planck ERCSC data also suggest that the original 
accelerated electron energy spectrum could be much harder than commonly thought, with power-law index around 1.5 instead of the canonical 
2.5. The implications of this are discussed for the acceleration mechanisms effective in blazar shock. Furthermore in many cases the Planck data 
indicate that gamma-ray emission must originate in the same shocks that produce the radio emission. 

Key words, galaxies: active - BL Lacertae objects: general - quasars: general - radiation mechanisms: non-thermal - radio continuum: galaxies 
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1. Introduction 

This paper is part of the first series of publications based 
on measurements made with the PlancJu satellite. Planck 



1 Planck (http://www.esa.int/planck) is a project of the European 

* Corresponding author: A. Lahteenmaki, alien@kurp.hut.fi Space Agency (ESA) with instruments provided by two scientific con- 
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dTauber et al.1 l201Qt iPlanck Collaboration! l20TTah is the third- 
generation space mission to measure the anisotropy of the cos- 
mic microwave background (CMB). It observes the sky in nine 
frequency bands covering 30-857 GHz with high sensitivity 
and angular res olution from 31/ to 5' . The Low Frequency 
Instrument (LFI: iMandolesi et al.l [20ld: iBersanelli et al.1 1201(1 
iMennella et al.ll201 ll) covers the 30. 44, and 70 GHz bands with 
amplifiers cooled to 20 K. The High Frequency Inst rument (HFI; 
lLamarre et aD 120101: IPlanck HFI Core Team! [20 1 1 ah covers the 
100, 143, 217, 353, 545, and 857 GHz bands with bolome- 
ters cooled to 0.1 K Polarization is measured in all but the 
highest two bands dLeahv et al.1 120101: iRosset et al] l201Qh . A 
combination of radiative cooling and three mechanical cool- 
ers produces the temperatures ne eded for the detectors and op- 
tics (Pla nck Collaborationll201 lbl) . Two data processing centres 
(DPCs) check and calibrate the data and make maps of the sky 
(IPlanck HFI Core Team! l2011ri IZacchei et all 201 1). Planck's 
sensitivity, angular resolution, and frequency coverage make it a 
powerful instrument for Galactic and extragalactic astrophysics 
as well as cosmology. 

The paper uses data from the Planck Early Release Co mpact 
Source Catalogue (ERCSC: IPlanck Collaboration! l2011cl) . The 
ERCSC provides positions and flux densities of compact sources 
found in each of the nine Planck frequency maps. The flux densi- 
ties are calculated using aperture photometry, with careful mod- 
elling of Planck's elliptical beams. The ERCSC includes data 
from the first all-sky survey, taken between 13 August 2009 and 
6 June 2010. This unique dataset offers the first glimpse of the 
previously unmapped millimetre and sub-millimetre sky. It is 
used here to create spectral energy distributions of 104 radio- 
bright, northern active galactic nuclei (AGN), with the most 
complete coverage in the radio to sub-millimetre frequencies to 
date. 

Radio-loud AGN host jets of relativistic matter emanat- 
ing symmetrically from the core. These jets produce copious 
amounts of non-thermal radiation, which dominates the spec- 
tral energy distributions (SEDs) of such sources compared with 
any thermal emission from the nucleus, i.e., the accretion disk. 
The SEDs typically consist of two broad-band bumps, the one 
at lower frequencies attributed to synchrotron radiation, and the 
other at higher frequencies attributed to inverse Compton (IC) 
radiation. The peak frequencies of the two bumps vary from 
one object to another. The peak of the synchrotron component 
can be between the infrared and X-ray domains, and the IC 
peak can range from MeV to GeV energies. The sequence in 
the peak frequencies of the emitted energy, and the factors that 
create it, have been a hot topic in blazar research for more 
than a decade (e.g.. iFossati et al.l 119971: IGhisellini et al.l 119981: 
Padovanil 120071: IG hisellini & Tavecchio 2008; Niep pola et atl 



20081: IGhisellini & Tavecchi o 2008; Sambruna et al. 20l3). 

Research on blazar SEDs has concentrated on two ap- 
proaches. The first is fitting the SED with a simple function, 
usually a second or third order polynomial, to obtain the pivotal 
peak frequencies and luminosities of the radiation components in 
a straightforward mann er. This approach is typically used when 
studying large samples (Fossati et al. 1998; Nieppola et al. 2006; 
ISambruna etaDl2006l: Nieppol a~al.l 120081: 1 Abdo et al.ll2010ah . 
The second approach is detailed modelling of the SED, starting 
with the definition of initial parameters such as electron energy, 
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magnetic field intensity, and Doppler factor. This method is more 
time-consuming and is used mostly for individual sources (e.g., 
lAcciari et al.ll2010l:ICollmar et al.ll2010l) . The standard model is 
a leptonic, homogeneous, one-zone model, where the emission 
originates in a si ngle componen t (for a review of the blazar emis- 
sion models, see Bottcher 2010). One-zone models are useful as 
first-order approximations, but, in reality, AGN jets are rarely, if 
ever, dominated by a single source of radiation. The material in 
the jets flows through shocks in the jet , which locally enhance 
the radiation (the "s hock-in- jet" model: iMarscher & Gear|[l985l: 
IValtaoja et al.ll 19921) . There can be several of these shocks in the 
jet simultaneously, and adding these to the emission of the quies- 
cent jet, we have several radiation components. Therefore, ide- 
ally, the SEDs and radio spectra should be modelled with more 
than one component. Such modelling is also necessary for the 
proper identification of the high frequency (IC) emission sites. 

In this paper we use the nearly complete SEDs and well- 
covered radio spectra provided by Planck to look for signs of 
these multiple components contributing to the total radiation that 
we observe. The paper is structured as follows. In Sects. [2] and [3] 
we introduce our sample and summarize the multifrequency data 
used in our study. The spectral energy distributions and radio 
spectra are discussed in Sects. |4] and [3 and their modelling is 
described in Sect. [6] In Sect. |7] we discuss the implications of our 
results for understanding the acceleration mechanisms in blazar 
jets, and in Sect. [8] we summarize our conclusions. Throughout 
the paper we adopt the sign convention for spectral index, a: 
S v oc v a . The errors of numerical values marked with a plus- 
minus sign correspond to one standard deviation. 



2. The sample 

The complete sample presented in this paper consists of 104 
northern and equatorial radio-loud AGN. It includes all AGN 
with declination > -10° that have a measured average radio flux 
density at 37 GHz exceeding 1 Jy. Most of the sample sources 
have been monitored at Metsahovi Radio Observatory for many 
years, and the brightest sources have been observed for up to 
30 years. The sample can be divided into subclasses as fol- 
lows: 33 high-polarization quasars (HPQs), 21 low-polarization 
quasars (LPQs), 21 BL Lacertae objects (BLOs), 19 normal 
quasars (QSOs), 9 radio galaxies (GALs), and one unclassi- 
fied source (J184915+67064). (See, for example. iHovatta et al.1 
d2008ll2009h for additional information on the classification.) By 
high-polarized quasars we mean objects which have a measured 
optical polarization > 3% at some epoch, while low-polarized 
quasars have a polarization < 3%. Normal quasars have no po- 
larization measured, so they could be either HPQs or LPQs. 
Radio galaxies are non-quasar AGN. The full sample is listed 
in Table [T] Columns 1 and 2 give the name and J2000 name for 
the source, and for some sources an alternative name is given in 
Col. 3. The coordinates of the sources are given in Cols. 4 and 
5. The start dates of the Planck scans are listed in Cols. 6 and 7. 
The average 37 GHz flux density from Metsahovi observations 
is given in Col. 8. For Col. 9 onwards, see Sect.ffl 



3. Multifrequency data 

The core of our data set is the Planck ERCSC. The con- 
struction and contents of the catalogue are described in 
Plan ck Collaboration! (1201 id) . For most sources, the ERCSC 
flux density values are averages of two scans, separated by about 
six months. To enable extensive multifrequency studies with a 
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simultaneous data set, the Planck Extragalactic Working Group 
has coordinated a programme in which ground-based and space- 
borne telescopes observe the sources in unison with Planck. In 
this paper we present SEDs based on the averaged ERCSC flux 
densities and on multifrequency supporting observations taken 
within two weeks of the Planck scans. In a later paper we will 
present single-epoch SEDs, constructed from all the available 
Planck and supplementary simultaneous multifrequency data. 
This collaborative programme includes 12 observatories around 
the world (Table |2]). Archival data have been obtained from 
the literature and from the the search tool at the ASI (Agenzia 
Spaziale Italiana) Science Data Center (ASDC) web page0. 

3.1. Radio and submillimetre data 

Centimetre-band observations were obtained with the University 
of Michigan Radio Astronomy Observatory's (UMRAO) 26- 
m prime focus paraboloid equipped with radiometers operating 
at central frequencies of 4.8, 8.0, and 14.5 GHz. Observations 
at all three frequencies employed rotating polarimeter systems 
permitting both total flux density and linear polarization to be 
measured. A typical measurement consisted of 8-16 individ- 
ual measurements over a 20-40 min period. Frequent drift scans 
were made across stronger sources to verify the telescope point- 
ing correction curves, and observations of programme sources 
were intermixed with observations of a grid of calibrator sources 
to correct for temporal changes in the antenna aperture effi- 
ciency. The flux sc ale was based on observations of Cassiopeia A 
(iBaars et al.lll977l) . De tails of the calib ration and analysis tech- 
niques are described in lAUer et al.l (Il985l) . 

Six-frequency broadband radio spectra were obtained with 
the RATAN-600 radio telescope in transit mode by obs erving si- 
multaneously at 1.1, 2.3, 4.8, 7 .7, 11.2. and 21.7 GHz (iPariiskiil 
1993; Berlin & Friedman 1996). The parameters of the receivers 
are listed in Tabled where v c is the central frequency, Ay is the 
bandwidth, AT is the sensitivity of the radiometer over a 1 s in- 
tegration, r p hys is the physical temperature of the front-end am- 
plifier, and r sys is the noise temperature of the whole system at 
the given frequency. Data were reduced using the RATAN stan- 
dard software FADPS (Flexi ble Astronomical D ata Processing 
System) reduction package ( Verkhodanov 1997). The flux den- 
sity measurement pro cedure at RATAN-600 is described by 
lAliakberov etaD(ll985h . 

The 37 GHz observations were made with the 13.7-m 
Metsahovi radio telescope using a 1 GHz bandwidth, dual-beam 
receiver centred at 36.8 GHz. The observations were ON-ON 
observations, alternating the source and the sky in each feed 
horn. The integration time used to obtain each flux density data 
point typically ranged from 1200 to 1400 s. The detection limit 
of the telescope at 37 GHz is ~ 0.2 Jy under optimal conditions. 
Data points with a signal-to-noise ratio less than four were han- 
dled as non-detections. The flux density scale was set by obser- 
vations of DR21. Sources NGC7027, 3C274, and 3C84 were 
used as secondary calibrators. A d etailed description of th e data 
reduction and analysis is given in iTerasranta et al.l (Il998l) . The 
error estimate in the flux density includes the contribution from 
the measurement rms and the uncertainty of the absolute calibra- 
tion. 

Quasi- simultaneous centimetre/millimetre radio spectra for 
a large number of Planck blazars have been obtained within 
the framework of a Fermi-GST related monitoring pro- 
gramme of gamma-ray blazars (the F-GAMMA programme, 

2 http://www.asdc.asi.it 



iFuhrmann etafl 120071: lAngelakis etaD 120081) . The frequency 
range spans 2.64 GHz to 142 GHz using the Effelsberg 100- 
m and IRAM 30-m telescopes. The Effelsberg measurements 
were conducted with the secondary focus heterodyne receivers at 
2.64, 4.85, 8.35, 10.45, 14.60, 23.05, 32.00, and 43.00 GHz. The 
observations were performed quasi- simultaneously with cross- 
scans, i.e., slewing over the source position in the azimuth and 
elevation directions, with an adaptive number of sub- scans for 
reach i ng the desired sensitivi ty (for details, see Fuhr mann et al.l 
2008; Angel akTs et al.l [2008b . Pointing offset correction, gain 
correction, atmospheric opacity correction, and sensitivity cor- 
rection have been applied to the data. 

The Institut de Radioastronomie Millimetrique (IRAM) ob- 
servations were carried out with calibrated cross-scans using the 
EMIR horizontal and vertical polarisation receivers operating at 
86.2 and 142.3 GHz. The opacity-corrected intensities were con- 
verted to the standard temperature scale and corrected for small 
remaining pointing offsets and systematic gain-elevation effects. 
The conversion to the standard flux-density scale was done us- 
ing the instantaneous conversion factors derived from frequently 
observed primary calibrators (Mars, Uranus) and secondary cal- 
ibrators (W3(OH), K3-50A, NGC 7027). From this programme, 
radio spectra measured quasi- simultaneously with the Planck 
observations have been collected for a total of 37 Planck blazars 
during the period August 2009 to June 2010. 

Many of the sources in the sample were monitored at 
15 GHz using the 40-m telescope of the Owens Valley Radio 
Observatory (OVRO) as part of a larger monitoring programme 
(iRichards et al.l 1201 ll) . The flux density of each source was 
measured approximately twice weekly, with occasional gaps 
due to poor weather or instrumental problems. The telescope 
is equipped with a cooled receiver installed at prime focus, 
with two symmetric off-axis corrugated horn feeds sensitive to 
left circular polarization. The telescope and receiver combina- 
tion produces a pair of approximately Gaussian beams (157" 
FWHM), separated in azimuth by 12.'95. The receiver has a cen- 
tre frequency of 15.0 GHz, a 3.0GHz bandwidth, and a noise- 
equivalent reception bandwidth of 2.5 GHz. Measurements were 
made using a Dicke- switched dual-beam system, with a second 
level of switching in azimuth to alternate source and sky in each 
of the two horns. Calibration was referred to 3C 286, for which a 
flux density of 3.44 Jy at 15 GHz is assumed (iBaars et al.ll 19771) . 
D etails of the observatio ns, calibration, and analysis are given 
bv lRichardsetaD(l201lh . 

The Very Large Array (VLA) and (since spring 2010) the 
Expanded VLA also observed a subset of the sources as simulta- 
neously as possible. Most of the VLA and EVLA runs were brief 
1-2 hour blocks of time. In a one-hour block of time, in addition 
to flux calibrators and phase calibrators, typically 5-8 Planck 
sources were observed. In many cases, VLA flux density and 
phase calibrators were themselves of interest, since they were 
bright enough to be detected by Planck. For these bright sources, 
the integration times could be extremely short. At 4.86 GHz and 
8.46 GHz, each target was observed for approximately 30 s. The 
integration times were typically 100 s at 22.46 GHz and 120 s 
at 43.34 GHz. All VLA/EVLA flux density measurements were 
calibrated using standard values for one or both of the primary 
calibrator sources used by NRAO (3C 48 or 3C 286), and the uv 
data were flagged, calibrated and imaged using standard NRAO 
software, MPS or CASA. The VLA and EVLA were in different 
configurations at different times in the several months duration 
of the observations. As a consequence, the angular resolution 
changed. In addition, for a given configuration, the resolution 
was much finer at higher frequencies. For that reason, sources 
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Table 2. Optical and radio observatories that participated in the Planck multifrequency campaigns 



Radio 




Observatory 


Frequencies [GHz] 


APEX, Chile 


345 


ATCA, Australia 


4.5 - 40 


Effelsberg, Germany 


2.64 - 43 


IRAM Pico Veleta, Spain 


86, 142 


Medicina, Italy 


5,' 8.3 


Metsahovi, Finland 


37 


OVRO, USA 


15 


RATAN-600, Russia 


1.1,2.3,4.8, 7.7, 11.2,21.7 


UMRAO 


4.8, 8.0, 14.5 


VLA, USA 


5, 8, 22, 43 


Optical 




Observatory 


Band 


KVA, Spain 


R 


Xinglong, China 


i 



Table 3. Parameters for the RATAN-600 receivers. 



fc [GHz] 


Af [GHz] 


AT [mK] 


^phys [K] 


T sys [K] 


21.7 


2.5 


3.5 


15 


77 


11.2 


1.4 


3 


15 


65 


7.7 


1.0 


3 


15 


62 


4.8 


0.9 


2.2 


15 


39 


2.3 


0.4 


8 


310 


95 


1.1 


0.12 


15 


310 


105 



that showed signs of resolution in any VLA configuration or at 
any VLA frequency have been carefully flagged. 

The Planck-AICA Co-eval Observations (PACO) project 
dMassardi et alJ [20lTah consists of observations of a compila- 
tion of source s selected in the Austra lia Telescope 20 GHz sur- 
vey (AT20G; iMassardi et al.l l2011bl) taken with the Australia 
Telescope Compact Array (ATCA) in the frequency range 4.5- 
40 GHz. The observations were carried out at several epochs 
close in time with Planck satellite observations covering July 
2009 - August 2010. The PACO sample is a complete, flux- 
density limited, and spectrally- selected sample over the whole 
southern sky, with the exception of the region with Galactic lat- 
itude \b\ < 5°. 

The Simultaneou s Medicina Planck Experiment (SiMPlE; 
IProcopio et al.l l201lb uses the 32-m Medicina single-dish an- 
tenna at 5 and 8.3 GHz to observe the 263 sources of the 
New Extragalactic W MAP Point Source (NEWPS) sample 
dMassardi et al.l [2009b with S > 0°, and partially overlapping 
with the PACO observations for -10° < S < 0°. The project be- 
gan in June 2010. Because of the lack of simultaneity with the 
ERCSC, the data have been used here to match the PACO ob- 
servations in the overlapping region and to add information for 
S > 70°, a region so far poorly covered at 5 GHz. 

Twelve sources from our sample were observed in the 
submillimetre domain with the 12-m Atacama Pathfinder 
Experiment (APEX) in Chile. The observations were made us- 
ing the LABOCA bolometer array centred at 345 GHz. Data 
were taken at two epochs in 2009: September 3-4 2009 and 
November 12 2009. The data were reduced using the script pack- 
age minicrusJ^ version 30-Oct-2009, with Uranus used as cali- 
brator. 



3 http ://w w w. submm. caltech. edu/~ share/crush/ 



3.2. Optical data 

The optical observations were made with the 3 5 -cm KVA 
(Kungliga Vetenskapsakademien) telescope on La Palma, 
Canary islands. All observations were made through the 7?-band 
filter (AqS = 640 nm) using a Santa Barbara ST- 8 CCD camera 
with a gain factor of 2.3 e~/ADXJ and readout noise of 14 elec- 
trons. The binning of pixels by 2 x 2 pixels resulted in a plate 
scale of 07 98 per pixel. We obtained three to six exposures of 
180 s per target. The images were reduced in the standard way 
of subtracting the bias and dark frames and dividing by twilight 
flat-fields. 

The flux densities of the target and three to five stars in the 
target field were measured with aperture photometry and the 
magnitude difference between the target and a primary reference 
star in the same field was determined. Using differential mode 
makes the observations insensitive to variations in atmospheric 
transparency and accurate measurements can be obtained even 
in partially cloudy conditions. The 7?-band magnitude of the 
primary reference star was determined from observations made 
on photometric nig hts, using comparison stars in known blazar 
fields as calibrators dFiorucci & Tostilll99d:lFiorucci et al.ll 19981 : 
iRaiteri et alJll998b IVillata et alJll998b iNilsson et alJl2007h and 
taking into account the color term of the 7?-band filter employed. 
After the 7?-band magnitude of the primary reference star was de- 
termined, the object magnitudes were computed from the mag- 
nitude differences. At this phase we assumed V - R = 0.5 for the 
targets. Several stars in the field were used to check the quality 
of the photometry and stability of the primary reference. The un- 
certainties in the magnitudes include the contribution from both 
measurement and calibration errors. 

The monitoring at Xinglong Station, National Astronomical 
Observatories of China, was performed with a 60/90-cm f/3 
Schmidt telescope. The telescope is equipped with a 4096x4096 
E2V CCD, which has a pixel size of 12yum and a spatial reso- 



4 



Planck Collaboration: SEDs and radio spectra of northern AGN 



lution of T/3 per pixel. The observations were made with an i 
filter. Its central wavelength and passband width are 668.5 and 
51.4nm, respectively. The exposure times are mostly 120 s but 
can range from 60 to 180 s, depending on weather and moon 
phase. 

3.3. Swift 

The S wift Gamma-Ray-Burst (GRB) Explorer ( Ge hrels et al.l 
2004) is a multi-frequency space observatory devoted to the dis- 
covery and rapid follow up of GRBs. There are three instruments 
on board the space craft: the UV and Optical Teles cope (UVOT, 
Roming et ai1l2005l) . the X-Ray Telescope (XRT. iBurrows et al.1 
2005]) sensitive i n the 0.3 - 10.0 keV ba nd, and the Burst Alert 
Telescope (BAT, Barthelm y et al.l 12005) sensitive in the 15 - 
150keV band. Although the primary scientific goal of the satel- 
lite is the observation of GRBs, the wide frequency coverage is 
useful for AGN studies because it covers the region where the 
transition between synchrotron and inverse Compton emission 
usually occurs. 

When not engaged in GRB observations, Swift is available 
for target of opportunity (ToO) requests. The Swift team de- 
cided to devote an average of three ToO per week to the si- 
multaneous observations of Planck AGN. The simultaneity of 
Swift observations within two weeks to the Planck first or sec- 
ond scan, or both, is shown in Tabled marked sim_lst, sim_2nd, 
or sim_lst_2nd, respectively. 

3.3.1. UVOT 

Swift UVOT observations were carried out using the "filter of 
the day", i.e., one of the lenticular filters (V, B, U, UVWl, 
UVM2, and UVW2), unless otherwise specified in the ToO re- 
quest. Therefore images are not always available for all filters. 

The photometry analysis of all our sources was performed 
using the standard UVOT software distributed within the 
HEAsoft 6.8.0 package, and the calibration included in the latest 
release of the "Calibration Database" (CALDB). A specific pro- 
cedure has been developed to process all the ToO observations 
requested by the Planck project. Counts were extracted from an 
aperture of 5" radius for all filters and converted to flux densi- 
ties using the standard zero points (iPoole et al.ll2008l) . The flux 
densities were then de-reddened usi ng the appropriate val ues of 
E(B - V) for each source taken from Schlegel et al.l (Il998l) . with 
A^/E(B - V) ratios calculated for UVO T filters using the m ean 
galactic interstellar extinction curve from Fitzpatrick (1999). No 
variability was detected within single exposures in any filter. 

3.3.2. XRT 

The XRT is usually operated in "Auto State" mode which auto- 
matically adjusts the CCD read-out mode to the source bright- 
ness, in an at tempt to avoid pile-up (see iBurrows et al.l 12005: 
iHill et al.ll2004l for details of the XRT observing modes). As a 
consequence, part of the data sample was collected using the 
most sensitive Photon Counting (PC) mode while Windowed 
Timing (WT) mode was used for bright sources. 

The XRT data were processed w ith the XRTDAS software 
package (v. 2.5.1. ICapalbi et al.ll2005l) developed at the ASDC 
and distributed by the NASA High Energy Astrophysics Archive 
Research Center (HEASARC) within the HEASoft package (v. 
6.9). Event files were calibrated and cleaned with standard filter- 
ing criteria with the xrtpipeline task, using the latest calibration 



files available in the Swift CALDB. Events in the energy range 
0.3-10keV with grades 0-12 (PC mode) and 0-2 (WT mode) 
were used for the analysis. 

Events for the spectral analysis were selected within a cir- 
cle of a 20-pixel (~ 47") radius, centred on the source posi- 
tion, which enclo ses about 90% of the point spread function 
(PSF) at 1.5 keV dMoretti et al.ll2005l) . For PC mode data, when 
the source count rate is above ~ 0.5 counts/s data are signifi- 
cantly affected by pile-up in the inner part of the PSF. For such 
cases, after comparin g the observed PSF p rofile with the analyt- 
ical model derived by Mor etti etaD (120051) . we removed pile-up 
effects by excluding events detected within a certain inner radius 
around the source position, and used an outer radius of 30 pix- 
els. The value of the inner radius was evaluated individually for 
each observation affected by pile-up, depending on the observed 
source count rate. 

Ancillary response files were generated with the xrtmkarf 
task applying corrections for the PSF losses and CCD defects. 
Source spectra were binned to ensure a minimum of 20 counts 
per bin to allow the x 2 minimization fitting technique to be 
used. We fitted the spectra adopting an absorbed power-law 
model with photon index T x . When deviations from a single 
power-law model were found, we adopted a log-parabo lic law 
of the form F(E) = KE~ {a+blo ^ E)) dMassaro et al.ll2004 which 
has been shown to fit the X-ray spectrum of blazars well (e.g., 
iGiommi et al . 1120051 iTramacere et al.ll2009l) . This spectral model 
is described by only two parameters: a, the photon index at 
1 keV, and b, the curvature of the parabola. For both models 
the amount of hydrogen-equivalent column den sity (Nu) was 
fixed to the Galactic value along the line of sight ( Kal berla et al.l 
120051) . 

3.4. Fermi-LAT Observations and Data Analysis 

The Large Area Telescope (LAT) onboard Fermi is an electron- 
positron pair-conversion telescope sensitive to gamma-rays of 
energies from 20 MeV to above 300 GeV The Fermi-LAI con- 
sists of a high-resolution silicon microstrip tracker, a Csl ho- 
doscopic electromagnetic calorimeter, and an anticoincidence 
detector for charged particles background identification. A full 
d escription of the instr ument and its performance can be found 
in lAtwood et al.l (120091) . The large field of view (~2.4 sr) allows 
the LAT to observe the full sky in survey mode every 3 hours. 
The LAT PSF strongly depends on both the energy and the con- 
version point in the tracker, but less on the incidence angle. For 
1 GeV normal incidence conversions in the upper section of the 
tracker, the PSF 68% containment radius is 0.6° . 

The Fermi-LNX data considered for this analysis cover the 
period from 4 August 2008 to 4 November 2010. They have been 
analyzed using the standard Fermi-LAI ScienceTools software 
packag^3 (version v9rl6) and selecting for each source only 
photons above 100 MeV, be longing to the diffuse class (Pass6 
V3 IRF, lAtwood et al.ll2009b which have the lowest background 
contamination. For each source, we selected only photons within 
a 15° region of interest (Rol) centred around the source itself. 
In order to avoid background contamination from the bright 
Earth limb, time intervals where the Earth entered the LAT Field 
of View (FoV) were excluded from the data sample. In addi- 
tion, we have excluded observations in which the source was 
viewed at zenith angles larger than 105°, where Earth's albedo 
gamma-rays increase the background contamination. The data 
were analyzed with a binned maximum likelihood technique 

4 http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
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(iMattox et al.| [T996) using the analysis software (gtlike) devel- 
oped by the LAT teairQ. A model accounting for the diffuse 
emission as well as the nearby gamma-ray sources is included 
in the fit. 

The diffuse foreground, including Galactic interstellar emis- 
sion, extragalactic gamma-ray emission, and residual cos- 
mic ray (CR) background, has been modelled using the 
model^j gllJem.v02 for the Galactic diffuse emission and 
isotropic Jem_v02 for the extra-galactic isotropic emission. Each 
source has been fit with a power law function 

dN N(y+1)E7 



where both the normalization factor N and the photon index y 
are left free in the model fit. The model also includes all the 
sources wi thin a 20° Rol inc luded in the Fermi-LAI one year 
catalogue (lAbdo et al.ll2010bl) . modelled using power law func- 
tions. If a source include d in the model is a pulsar belonging to 
the Fermi pulsar catalog (lAbdo et al.ll2010dl) . we have modelled 
the source using a power-law with exponential cut-off and the 
spectral parameters from the pulsar catalogue. For the evaluation 
of the gamma-ray SEDs, the whole energy range from 100 MeV 
to 300 GeV is divided into two equal logarithmically spaced bins 
per decade. In each energy bin the standard gtlike binned anal- 
ysis has been applied assuming for all the point sources in the 
model a power law spectrum with photon index fixed to -2.0. 
Assuming that in each energy bin the spectral shape can be ap- 
proximated by a power law, the flux density of the source in all 
selected energy bins was evaluated, requiring in each energy bin 
a test statistic (TSfl greater than 10. If the TS is lower than 10, an 
upper limit is evaluated in that energy bin. Only statistical errors 
for the fit parameters are shown in the plots. Systematic effects 
are mainly based on uncertainties in the LAT effective area de- 
rived from the on-orbit estimations, and are <5% near 1 GeV, 
10% below 0.10 GeV, and 20% above 10 GeV. 

The LAT gamma-ray spectra of all AGN detected by Fermi 
are studied in lAbdo e t al. ( 2010c), based on 1 1 months of Fermi- 
LAT data. For this paper, we derived the gamma-ray SEDs 
of the 104 sources in the sample in three time intervals, pre- 
sented in Table H In the first interval the Fermi-LAI observa- 
tions are simultaneous to the Planck first or second scan, or both, 
within two weeks (marked sim_lst, sim_2nd, or sim_lst_2nd in 
Table |4]). The Fermi-LAI data have in this case been integrated 
over two weeks. In the second interval the gamma-ray data are 
quasi- simultaneous to the Planck first or second scan, or both, 
having been integrated over two months (2M_lst, 2M_2nd, or 
2M_lst_2nd in Table [4]). In the third interval are sources for 
which Fermi-LAI data have been averaged over 27 months due 
to their faintness (27m in Table [4]). 

Note that eight sources (0306+102, 0355+508, 0804+499, 
0945+408, Jl 130+3815, 1413+135, 1928+738, 2005+403) that 
appear as significant gamma-ray detections in the 27 month 
dataset presented in this work, h ave not been inclu ded in any 
Fermi catalogue published so far (lAbdo et al.ll2010bl) . They can 
be considered new gamma-ray emitting sources, taking into ac- 



5 http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
Cicerone_Likelihood 

6 http ://f ermi . gsf c . nasa. go v/s sc/data/acces s/lat/ 
BackgroundModels.html 

7 The test statistic (TS) is defined as TS = -21n(L /Li) with L the 
likelihood of the null-hypothesis model as compared to the likelihood 
of a competitive model, L\ 



count the possibility that their association with the radio source 
is purely spatial. 

4. Spectral energy distributions 

The SEDs for the whole sample are shown in Figs. [T8l - 11211 (in 
panels on the left). The Planck measurements and the simulta- 
neous auxiliary data are indicated by red circles. The grey points 
represent archival data obtained from the literature and using the 
ASDC search tool. The measurements were regarded as simul- 
taneous if they were taken within two weeks (for radio frequen- 
cies) or five days (for optical and higher energies) of the Planck 
measurement. Previous studies of the radi o variability of blazars 
dHovatta et al.ll2008l : iNieppola et"aDl2009l) have shown that the 
two-week simultaneity limit is appropriate, as large-scale radio 
flux-density changes on that timescale are quite rare. The SEDs 
were fitted with second and third degree polynomials, namely 

log vF v = ci(log v) 2 + c 2 (log v) + c 3 , (2) 

and 

log vF y = c 4 (log y) 3 + c 5 (log y) 2 + c 6 (log v) + c 7 , (3) 

where ci (i = 1 ... 7) are fit parameters. We calculated the peak 
frequencies of the components from the third degree fit by de- 
fault and used the second degree fit only if it seemed by visual in- 
spection to be more accurate. Typically, in these cases, the third 
degree fit would not yield a sensible result at all. In the SEDs in 
Figs. [T8l - 11211 third degree fits are marked with a dashed line, 
second degree fits with a dotted line. 

The fits were divided into two classes, A and B, according 
to the fit quality, A being superior. We emphasize that this divi- 
sion is the subjective view of the authors, and is only meant to 
illustrate the varying quality of the SED fits. The classes have 
not been used in calculating the correlation coefficients (see be- 
low). We were able to obtain meaningful synchrotron fits for 60 
sources (15 class A fits), and IC fits for 30 sources (10 class 
A fits). For 21 sources we were able to fit both components. 
The synchrotron and IC peak frequencies are listed in Table [T] 
Columns 9 and 10 give the peak frequency and the fit quality for 
the synchrotron component, and Cols. 11 and 12 give the peak 
frequency and the fit quality for the inverse Compton compo- 
nent. 

Figure[T] shows the distributions of the logarithms of the peak 
frequencies. The synchrotron peak frequencies, logyjj, are typi- 
cally very low, as we are dealing with bright radio sources. They 
range over two orders of magnitude, from 12.2 to 14.3, the av- 
erage of the distribution being 13.2 + 0.5. The source with the 
highest yjj in the near-ultraviolet domain is 0716+714. The IC 
peak frequencies, log y p c , range over three orders of magnitude, 
from MeV to the high GeV gamma-ray region. The distribution 
average is 22.0 + 0.9. The source with the highest IC peak fre- 
quency is 1156+295, with logy* = 23.5. The source with the 
lowest IC peak frequency, with log v* c = 20.5, is 0836+710. 

We have plotted the interdependence of the component peaks 
in Fig.O The two peak frequencies seem to have no significant 
correlation, which is confirmed by the Spearman rank correla- 
tion test (p = 0.301, P = 0.092, both class fits included). There 
is, however, a tendency for the sources with high yjj to have a 
high yjp. In the synchrotron- self-Compton (SSC) scenario, the 
separation of the component peaks depends on many factors, 
such as the electron Lorentz factors, magnetic field strengths, 
and particle densities. Therefore the lack of correlation as such 
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Fig. 1. The distributions of the synchrotron (top panel) and IC 
(bottom panel) peak frequencies. 



does not allow any strong statements to be made about the origin 
of the IC radiation (SSC, or external Compton, EC). 

We also compared our peak frequencies with those of 
Abdo et al.1 (l2010ak . These comparisons are illustrated in Fig. [3 
where the one-to-one correspondence is shown with a dashed 
line. Both correlations are significant according to the Spearman 
test (p = 0.629, P = 0.003 for the synchrotron peaks and 
p = 0.660, P = 0.010 for the IC ones). 



5. Radio spectra 

In addition to the spectral energy distributions, we have plotted 
the standard radio spectra for the sample sources. The aim was 
to look for signs of multiple components contributing to the syn- 
chrotron emission. The spectra are shown in Figs. [T814121I (in 
panels on the right). Red circles indicate low frequency (LF; < 
70 GHz) data simultaneous to Planck, red stars indicate ERCSC 
LFI data, blue circles indicate high frequency (HF; > 70 GHz) 
data simultaneous to Planck, and blue stars indicate ERCSC HFI 
data. As can be seen, the Planck data and the supporting multi- 
frequency data provide an excellent opportunity to find the sig- 
natures of possible multiple components or anomalous spectral 
shapes. A statistical study of extragalactic radio sources in the 
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ERCSC is presented in lPlanck Collaboration! (1201 lib , and SEDs 
of extreme radio sources, such as Gigahertz- Peaked Spectrum 
(GPS) sources in lPlanck Collaboration! (1201 lib . 

In Figs. |4] and [5] we have plotted the distributions of the LF 
and HF spectral indices using only the Planck data for the latter. 
When plotting the HF a values, we have only taken into account 
those sources with three or more data points. This leaves us with 
84 sources out of the 104. The sources with only one Planck 
scan are shown hatched in Figs. |4] and Figure [6] shows c#lf 
versus auF- As expected, the LF indices are fairly flat, with an 
average of -0.06. Their distribution is narrow, with 91% of the 
indices being in the range c*lf = -0.5-0.5. There are a couple 
of sources with remarkably steep LF spectra, namely 0552+398 
and 2021+614, both having Qf LF < -0.8. At the other end we 
have 0007+106 which has an inverted spectrum with c*lf = 0.86. 
Also 1228+126 has a steep spectrum, although the fit value may 
be exaggerated in this case (see Fig. [75]). For clarity, this source 
has been omitted from Figs. |4] and [6] 

The HF indices are concentrated around -0.5, having an av- 
erage of -0.56. Rather simila r distributions for LF and H F spec- 
tral indices were found in iPlanck Collaboration! (1201 lil) . espe- 
cially in their Fig. 7. However, a large part of the auF distribu- 
tion is in the flat domain with auF > -0.5. Extreme examples are 
1413+135 (a U F = 0.02) and 0954+658 (a HF = 0.34). It is unex- 
pected that the spectral indices of blazars at sub-millimetre and 
FIR frequencies should be this flat. Only 28 of the 84 sources 
with HF spectral fits shown in Figs. [T814121I have a < -0.7, 
the "canonical" value for optically thin spectra, corresponding 
to an electron energy index s of about 2.4. This has two possible 
explanations: either the total HF spectra are defined by several 
underlying components, or the energy spectrum of the electron 
population producing the radiation is much harder than generally 
assumed. Although some HF spectra show evidence for several 
sub-components or even an upturn at the highest frequencies, 
many others appear to be perfectly straight and therefore per- 
haps originate from a single optically thin component. At the 
highest Planck frequencies the lifetimes of radiating electrons 
are short, and one would in most cases expect to see the opti- 
cally thin spectra steepened by energy losses, with indices rang- 
ing from the original Qf t hin to Of t hin - 0.5 and beyond. This makes 
spectral indices > -0.7 even more remarkable. 

To rule out contamination by dust in the host galaxies, we 
have estimated the possible contribution of an extremely bright 
submillimetre galaxy. We also estimated the flux densities of 
2251 + 158 at several submillimetre bands in the case its lumi- 
nosity were similar to the submillimetre galaxy. The compari- 
son shows that if the dust component in 2251 + 158 were of the 
same brightness as in the extremely bright submillimetre galaxy, 
its emission would be completely swamped by the non-thermal 
emission from the jet. The LF and HF spectral indices have 
been calculated using ERCSC data which for most sources is 
an average of two Planck scans separated by about six months. 
However, the HF spectral index distribution for the sources that 
have been observed only once by Planck is quite similar to that 
of the whole sample. The implications of this result are discussed 
in more detail in Sect. H 

As Fig.[6]shows, BLOs, HPQs and LPQs have different ckhf 
distributions. In Table [5] we have listed the average spectral in- 
dices for both frequency ranges for all subgroups. The LF spec- 
tra of radio galaxies are very steep on average, but the value 
is greatly influenced by the uncertain fit of 1228+126. If this 
source is left out of the calculations, the average q^lf for galax- 
ies is 0.10 + 0.52. LPQs have the steepest spectra in the HF 
range. However, the standard deviations of the samples are sub- 



Table 5. Average spectral indices for AGN subclasses. The er- 
rors correspond to one standard deviation. 



Class 


number LF 


number HF 


average Qf LF 


average or H F 


BLO 


21 


20 


-0.01 + 0.25 


-0.34 + 0.26 


HPQ 


33 


28 


-0.06 + 0.22 


-0.58 + 0.27 


LPQ 


21 


18 


-0.08 + 0.30 


-0.80 + 0.27 


QSO 


19 


10 


-0.05 + 0.28 


-0.57 + 0.26 


GAL 
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-0.25 + 0.66 


-0.55 + 0.19 
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Fig. 4. Distribution of LF spectral indices for the whole sample 
(104 sources). The sources that have only been scanned once by 
Planck are shown hatched (46 sources). 

stantial, which can also be seen as the large scatter in Fig. [6] 
To quantify the possible differences between the AGN classes, 
we ran the Kruskal-Wallis ANOVA test. For the LF indices we 
could find no significant differences. However, in the case of the 
HF indices, the distributions of the subgroups differ significantly 
(P < 0.001). Multiple-comparisons test with t distribution tells 
us that LPQs differ from all other subgroups (P < 0.031). Also, 
BLO indices have a different distribution from all quasar sub- 
groups (P < 0.019). 

6. Modelling the SEDs and spectra of blazars 

Spectral energy distributions rapidly became one of the main 
tools for understanding blazar physics after the EGRET instru- 
ment aboard the Compton Gamma-ray Observatory (CGRO) 
satellite discovered strong gamma-ray emission from 3C279 
(Hartma n et al.ll 19921) and, subsequently, from a large number of 
other radio-bright AGN. It was rapidly accepted that the double- 
peaked overall shape of the radio-to-gamma SED was due to 
synchrotron and inverse Compton radiation, but beyond this the 
agreement ended. The data simply lacked the sensitivity and 
the time coverage to sufficiently constrain the mo dels. Even the 
most detailed and convincing effort for 3C 279 (lHartman et al.l 
200 lj) could be criti cized for unrealistic physical assumptions 
dLindfors et alJl2005b . 

With the coming of the Fermi and Swift satellites, ground- 
based TeV telescopes, and other satellite data, the inverse 
Compton high-energy SED of a large number of AGN can now 
be determined with hugely improved accuracy and time cov- 
erage. This has initiated a new era of AGN modelling (e.g., 
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Fig. 5. Distribution of HF spectral indices for the 84 sources that 
had three or more data points. The sources that have only been 
scanned once by Planck are shown hatched (40 sources). 
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Proceedings of the Work shop "Fermi meets Jansk y: AGN in 
Radio and Gamma Rays"; Fermi meets Janskv 2010). However, 
since the IC spectrum is created by the relativistic electrons scat- 
tering either the synchrotron photons (SSC) or external photons 
(EC), accurate knowledge of the shape of the synchrotron com- 
ponent is essential for all realistic modelling. Here the Planck 
data are invaluable, since they provide the first simultaneous 
multifrequency coverage of the crucial radio-to-IR part of the 
synchrotron spectrum. 

However, most of the ERCSC data are not suitable for de- 
tailed modelling, since they are averages of two spectra taken 
about six months apart. Considering the strong SED variabil- 
ity, accurate modelling can be done only with the final Planck 
datasets and their several individual scans. Here we limit our- 
selves to a short discussion of how Planck and supporting data 
can be used to improve our understanding of blazar physics. 

The basic picture of a blazar is well established: a black 
hole surrounded by an accretion disk and a broad-line emit- 



ting region (BLR), and a relativistic jet which produces the 
lower-frequency bump of the SED through synchrotron emis- 
sion. The jet is not totally stable; VLB I imaging shows how 
new components, identified as relativistic shocks in the jet, 
emerge from the radio core at intervals ranging from months 
to years. These growing and decaying shocks are also respon- 
sible for the total flux density variations in the radio-to-IR, 
and at least partly also in the optical regime. The important 
thing to note is that the synchrotron SED is rarely, if ever, the 
product of a single synchrotron-radiating component. The na- 
ture of the radio core is still being debated, as is the nature 
and the extent of the jet flow upstre am of the radio core, be- 
tween it and the bla ck hole (e.g., iMarscher & Jorstadl 1201 Ol: 
Marscher et al. 2 0To|). Shocked jet models (Marscher & Gear 
119851: iHughes et alJll989b IValtaoia et al1ll992h provide a good 
account of the multifrequency var iations in the total flu x density. 
In the particular case of 3C 279, iLindfors et al] (120061) showed 
that the optical variations can also be moderately well explained 
as the higher-frequency part of the shock emission. 

The origin of the IC component is a matter of intense de- 
bate. In the CGRO era, the favoured alternative was that the 
high-frequency radiation originates close to the black hole and 
the accretion disk, inside the BLR and well upstream of the ra- 
dio core. With a strong radiation field external to the jet, exter- 
nal Compton is then likely to be the dominant emission mech- 
anism. A minority view was that at least a part of the gamma- 
ray flaring comes from much farther out, at or downstream 
from the radio core. The evidence for this was that, statisti- 
cally speaking, strong gamma-ray flarin g tended to occur af- 
ter the beginning of a millime tre flare ( Valtaoja & Terasranta 
1 19961: lLahteenma ki & Valta oia 200 3]) and aft er the ejection of 
a new VLBI component dJorstad et al.ll200ll) . Occurring well 
outside the BLR, SSC seems in this case to be the only vi- 
able emission mechanism, but on the other hand simple SSC 
models have generally failed to explain the gamma-ray emis- 
sion, requiring unrealistically low magnetic field strengths in 
other than the extreme TeV BL Lac objects. With new Fermi 
data, the evidence for gamma-ray emi s sion from the radio jet 
has strengthened dKovalev et al] 120091: Pushkar ev et al] 12 1 ; 
Leon-Tavares et al. 2011; Niep pola et al.ll2010l:lTornikoski et al] 
l20ld:IValtaoja et al .120101) . but the source of the upscattered pho- 
tons remains a puzzle. 

Most of the attempts to model the high-energy bump have 
adopted a "reverse engineering" method, starting from the sec- 
ondary IC SED and not from the primary synchrotron SED. 
The original electron energy spectrum is considered to be rather 
freely adjustable, with parameters chosen to produce a good fit 
of the calculated IC spectrum to the observations within the as- 
sumed model. In most cases, the synchrotron spectrum is as- 
sumed to originate in a single homogeneous component, and of- 
ten the spectrum is adjusted only to the optical data, without 
any attempt to explain the radio-to-IR part of the SED. There 
is not much observational or theoretical evidence for the actual 
physical existence of such an IR-to-optical component: it is sim- 
ply assumed to be there, because it can produce the observed IC 
spectrum. 

To some extent such an approach is justifiable if the gamma- 
rays originate upstream of the radio core, and are therefore not 
directly connected to the observed radio jet and shocks, and their 
lower-frequency synchrotron emission. However, even in this 
case physically more realistic models must take into account also 
the radio-to-IR part, for several reasons. First, whatever the exact 
mechanism, the radio shocks also contribute to the IC spectrum, 
and their contribution must be accounted for in the model fits. 
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Assuming reasonable physical parameters, the shocks can quite 
often produce part or all o f the observed X-ray e mission through 
the SSC mechanism (e.g. jLindfors et al]|2006l) . Indeed, a stan- 
dard method for estimating the physical parameters of VLBI 
components i s to assume that the X-ray emission i s produced 
through SSC dGhisellini et al.lll993HMarscherll 19871) . Secondly, 
the components seen in the radio, especially in the Planck HFI 
regime, contribute to the optical synchrotron emission, and this 
contribution must be subtracted from the putative IR-to-optical 
component, which, as explained above, in many models gener- 
ates the whole IC spectrum. For example, in the present sam- 
ple, at least 0235+164, 0917+449, 0954+658, 2227-088, and 
especially 2251 + 158 show ongoing millimetre-to-submillimetre 
flaring, which leaves little room for any additional components 
between the highest Planck frequencies and the optical. This is 
the case in particular if the o ptical ly thin spectra are originally 
very flat, as we argue in Sects. l6Jl and 171 

On the other hand, if gamma-rays do come from the shocks 
at, or downstream from, the radio core, as increasing evidence 
indicates, the obvious first step is to model the shocks them- 
selves and their spectra as accurately as possible using Planck 
and other radio data, and only after that attempt to calculate the 
IC radiation they produce. 

Finally, it is clear that the more realistic multicomponent 
models, whatever their exact details, will give rather different 
physical parameter values than the single-component models. As 
an example, if the X-ray emission and the gamma-ray emission 
are not produced by a single component but by two (or even 
more) components, the derived electron energy spectral parame- 
ters will be quite different. 

6.1. Some examples of different spectra and their 
implications 

In this section we highlight five examples of different spectra of 
ERCSC sources which have been scanned only once by Planck 
and therefore have true snapshot spectra presented in Figs. [T8l - 
I121I As stated before, our aim at this stage is not to model these 
sources, but to draw attention to some features commonly seen 
in the spectral energy distributions, and to their implications for 
future multi-epoch modelling of the Planck sources. We have 
fitted the radio-to-optical spectra with one or several standard 
self-absorbed sychrotron components (based on a power-law as- 
sumption for the electron energy spectrum). These are not model 
fits from a numerical, physical code (except for Figs. [T6l and [TTt, 
but rather are meant to guide the eye and to illustrate the general 
shape of the radio-to-optical spectra. 

We also use support observations to address the behaviour 
of the five example sources. In Figs. I7UTT1 we have plotted the 
37 GHz long-term flux density curves for the five sources. The 
long-term data have been taken at Metsahovi Radio Observatory. 

Even without numerical modelling we can address two im- 
portant topics with these sources. The first one is the spectral 
flatness of the optically thin radiation possibly hinting at non- 
standard acceleration processes. The second one concerns the 
origin of the IC component. As Figs. [T2l - [T7l show there does not 
appear to be room for additional synchrotron components be- 
tween the Planck and the optical frequencies. The IC radiation 
must therefore originate in the same shocks that also produce the 
radio emission, in contradiction to most theoretical models. 

0234+285 (Fig. [T2l) . This source was in a very uneventful 
state during the Planck observing period (Fig. [7]) and shows a 
spectrum typical of sources in a quiescent stage. The radio spec- 
trum can be modelled with a single synchrotron component, hav- 
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Fig. 7. The 37 GHz lightcurve of 0234+285, measured at 
Metsahovi Radio Observatory, showing continuous variability. 
The period when Planck was observing the source is indicated 
by the hatched region in the lower panel, which is a blow-up of 
the last year. 
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Fig. 8. The 37 GHz lightcurve of 0235+164, details as in Fig.H 
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Fig. 9. The 37 GHz lightcurve of 1253-055 (3C 279), details as 
in Fig. 
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Fig. 10. The 37 GHz lightcurve of 1413+135, details as in Fig.[7] 
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Fig. 11. The 37 GHz lightcurve of 2251 + 158 (3C 454.3), details 
as in Fig. [7] 



ing arthin = -0.61, which apparently steepens by Aa = -0.5 
and smoothly joins the optical spectrum (see also Fig. [29]). This 
synchrotron component can be identified with the latest signifi- 
cant flare (new shock), which peaked about two years earlier at 
37 GHz (Fig. [7]). During low flux levels, in the absence of recent 
flares/shocks, the underlying, relatively stable jet may also con- 
tribute significantly to the synchrotron emission. Another similar 
example is the source 0007+106 (Fig.[T9b. Lacking mid-IR data, 
we cannot prove that there is no additional synchrotron compo- 
nent, between the Planck frequencies and the optical, that might 
be the main source of the IC flux (Fig. [29b. However, the most 
straightforward interpretation of the data is that the IC flux also 
comes from the synchrotron component outlined in Fig.[T2j from 
a rather old shock far downstream the jet, and far beyond the 
BLR. 

0235+164 (Fig. [13]). The bumpiness of the radio spectrum 
shows clear evidence of multiple components (shocks) con- 
tributing to the observed total flux density. We have outlined two 
possible synchrotron components, which together could explain 
most of the observed radio spectrum. Although 0235+164 was 
relatively quiescent during the Planck observation (Fig. [8]), there 
were at least four strong flares during the previous five years, and 
they all are likely to be still contributing significantly to the total 
radio flux density, at least around 10 10 Hz and below. Assuming 



reasonable physical parameters for these shocks, they can also 
produce significant amounts of X-rays through the SSC mecha- 
nism, which must be accounted for in a proper multicomponent 
modelling of the IC spectrum. As with 0234+285, the higher- 
turnover component with a±m = -0.4 can be joined smoothly to 
the optical spectrum assuming a steepening by about Aa = -0.5, 
leaving no room for additional mid-IR components. Fermi mea- 
sured only upper limits for the gamma-ray flux, which may be 
due to the fact that even the most recent shock, peaking around 
100 GHz, is already about two years old (Fig. [8]). 

1253-055 (3C279). (Fig. O. Again, the highest-peaking 
radio component with a±m = -0.6 can be made to smoothly 
join the steep optical spectrum. From Fig.[l4]it is rather obvious 
that there can be no additional mid-IR components which might 
provide a significant contribution to the IC spectrum. If we as- 
sume that there is another synchrotron component turning over 
somewhere above the Planck frequencies and joining smoothly 
to the optical spectrum, Fig. [14] shows that its flux density can- 
not exceed the dashed line at any IR frequency between 10 12 and 
10 14 5 Hz. This is the case even if we make the assumption that 
the spectrum of the component seen at the highest Planck fre- 
quencies steepens very rapidly above 857 GHz. Consequently, 
the Fermi gamma-rays must originate in the latest shock, which 
started a few months before the Planck observation (Fig. [9]). 

1413+135 (Fig. [15]). This source exhibits a totally flat spec- 
trum, with both the LF and the HF spectral indices approxi- 
mately 0. Possible explanations are: an underlying electron en- 
ergy index s ~ 1, which is rather unlikely; the superposition 
of many self-absorbed synchrotron components (shocks of var- 
ious ages) with different turnover frequencies; or a continuous 
jet dominating the flux and having an integrated spectral index 
close to zero over a wide frequency range. As Fig. [TO] shows, the 
source has been rather inactive for over a decade, with no ma- 
jor ourbursts. The HFI spectrum may indicate an upturn, hinting 
at the presence of an IR component. We show a possible spec- 
tral decomposition, including an IR-peaking component. Further 
Planck observations are needed to ascertain the reality and the 
nature of this component, which may be either a synchrotron 
flare seen at an early developmental stage or, possibly, an in- 
frared dust component. The correct alternative can be identified 
through variability information. 

2251+158 (3C 454.3). (Figs. [TO] and [[7]). This source ex- 
hibits the most spectacular case of ongoing IR flaring, show- 
ing a very strong synchrotron component with a self- absorption 
turnover at about 80 GHz and a very flat (a = -0.2) optically 
thin spectrum. As Fig. Qj] shows, a series of spectacular flar- 
ing events has occurred in 2251 + 158 during the last five years, 
with Planck observing the source during the early stages of the 
strongest outburst ever seen in the source. In Fig. [TO] we show 
a model fit to the radio-to-optical data, using a preliminary ver- 
sion of the code that we are developing (Tammi et al., in prepa- 
ration). A close-up of the fit in the radio regime is shown in 
Fig. [T7] As Fig. 11191 also shows, it is impossible to squeeze in 
another component between the strong radio-IR flare and the 
optical regime. The gamma-ray emission must come from the 
growing shock that we see both in the total flux density monitor- 
ing curve (Fig.fTTl) and in the SED (Fig. II 19b . Most of the theo- 
retical model scenarios assume that gamma-rays originate close 
to the black hole and the accretion disk, well within the BLR and 
far upstream of the radio core. This is not a viable alternative for 
225 1 + 158, in which gamma-rays must come from the new shock 
component, parsecs downstream from the radio core. As Fig.fTTI 
shows, at the time of the Planck observation the shock had al- 
ready been growing for about half a year, which, changed into 
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Fig. 12. A possible single-component spectrum for 0234+285 to Fig. 14. The radio spectrum of 1253-055 shows that there is no 
illustrate the general shape. Red circles, LF data simultaneous to room for another component between radio and optical regimes. 



Planck: blue circles, HF data simultaneous to Planck. 
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Fig. 13. A possible multiple-component spectrum for 0235+164 Fi §- 15. The radio spectrum of 1413+135 is very flat and can be 



to illustrate the general shape. Symbols as in Fig. [121 



fitted with multiple components. Symbols as in Fig.[ 



the source frame, translates into a distance L = D 2 c(l +z) _1 A^ b s 
(where D is the Doppler boosting factor, c the speed of light, 
z the redshift, and Ar b s the time elapsed from the onset of the 
flare), of at least 10 pc down the jet, assuming a m odest value of 
D ~ 1 0, or 90 pc using the value of D = 33.2 from lHovatta et al.l 
(2009). This is the distance that the shock has moved down- 
stream after emerging from the radio core. 



7. The shape of the synchrotron spectrum: 
nonstandard acceleration mechanisms? 

If the electron energy distribution is assumed to be a simple 
power-law with the form N(E) = KE~ S , the synchrotron spec- 
trum of a single homogeneous source has an optically thick low- 
frequency part with a spectral index a = +2.5, a turnover, and an 
optically thin part with a spectral index a related to the electron 
energy index by a = -(s - l)/2. Energy losses eventually cause 
the spectrum to steepen to a - 0.5. 



1 cr u i 




v [Hz] 

Fig. 16. The radio-optical SED of 2251 + 158 (3C 454.3) fitted 
with one component and synchrotron emission from the jet 
(dashed line) using a model by Tammi et al. (in preparation). 
The radio spectrum is shown in Fig.fTTl 
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Fig. 17. A close-up of the radio spectrum and model fit of 
2251 + 158 (3C 454.3), details as in Fig.[T5| 



Since the earliest times of AGN research, the canonical value 
for a has been assumed to be about -0.75, corresponding to s = 
2.5. Such an electron energy index emerges naturally from the 
simplest first-order Fermi acceleration mechanism, and blazars 
also tend to show slopes around this value at the high-frequency 
(10 to 100 GHz) end of the radio spectrum, where the sources 
are typically assumed to be optically thin. However, Planck data 
now challenge this commonly accepted view. As Fig. [4] shows, 
the LF spectral slopes are, in the main, quite flat, and sometimes 
even rising. It is likely that in most cases this results from the 
superposition of several distinct synchrotron components, since 
the spectra are generally not smooth, and we also know from 
VLBI observations that there are nearly always several compo- 
nents (shocks) present in any source (see also Sect.0). The HF 
spectra (Fig. O, on the other hand, are often straight within the 
Planck flux density error bars, and give a definite impression 
of optically thin radiation. In some cases, such as 2251 + 158, 
the spectrum even shows a clear self-absorption turnover at high 
frequencies, followed by a flat optically thin part with a = -0.2 

(Figs.H5|andQ3. 

The flatness of the HF spectra, shown in Fig. [5] is remark- 
able. Only a third of the sources have a steeper than -0.7, while 
15 have a > -0.3. While the present data are not sufficient to ex- 
clude the possibility that multiple components produce the flat- 
ness, and some sources such as 0954+658 do show evidence o f 
an additional IR component (see also, e.g.. iRaiteri et al.|[l999l) . 
we consider such an explanation unlikely for the whole sample. 
The errors in auF cannot explain the very flat spectral indices, 
since they generally are smaller than Aa « 0.1. 

In Table [6] we list the 10 sources that fullfill the following 
criteria. They have been scanned only once by Planck and they 
have auF > -0.5 + Aa, where Aa is the HF spectral fit error 
shown in Figs. [T8l- 11211 These form the most extreme end of 
the QfRF distribution shown in Fig. \5\ Three of them may show 
indications of multiple HF components in the spectra, but for the 
seven other sources the HF spectra appear to be straight within 
the Planck errors. If this is indeed the case, they cannot be ex- 
plained with standard acceleration mechanisms with s > 2. 

The rise of a submillimetre flare is quite fast (see, for exam- 
ple, the JCMT monitoring data of [Stevens et alJll994h : the self- 
absorption turnover of a fresh shock component passes through 
the Planck HFI frequencies rather rapidly, and one therefore 
sees it only rarely in the snapshot high frequency spectra. (In 
the present data only 2251 + 158 is clearly such a case, al- 



Table 6. A list of the 10 sources scanned only once by Planck 
and with auF > -0.5 +Aa. S = straight HF spectrum; B = bumpy 
HF spectrum, possibly indicating multiple components. 



Source 


Qf H F 


Spectral shape 


0003-066 


-0.38 + 0.09 


S 


0059+581 


-0.19 + 0.09 


s 


0235+164 


-0.40 + 0.03 


s 


0430+052 


-0.34 + 0.05 


s 


1308+326 


-0.41 + 0.08 


B? 


1413+135 


0.02 + 0.12 


B 


1418+546 


-0.16 + 0.04 


S 


1652+398 


-0.21 + 0.04 


B 


1823+568 


-0.31+0.13 


S 


2251+158 


-0.20 + 0.02 


s 



though 0235+164 and 1652+398 might be others.) Instead, in 
most cases we expect to see optically thin submillimetre spec- 
tra which have already experienced energy losses. Therefore the 
observed high-frequency spectral-index distribution of a sam- 
ple of sources should show a range from a±m = ~(s - l)/2 to 
<*thin -0.5 and beyond, as the exponential steepening sets in, with 
spectral indices steeper than the initial a±m dominating. The ob- 
served distribution in Fig. [5] is thus incompatible with an initial 
Qffhin around -0.7, as we should then see a few values around 
a = -0.7 and a distribution peaking towards -1.2, reflecting the 
Aa = -0.5 steepening due to energy losses. 

Instead, the spectral index distribution is compatible with an 
electron index s « 1.5, resulting in a distribution of a±m hav- 
ing smallest values around -0.2 to -0.3 and a maximum around 
-0.7. Such a suggestion of a very har d original electron spec- 
trum is not entirely new; IValtaoja et a l. (1988) stu died the shock 
spectr a and concluded that a is around -0.2, and Hu ghes et al.l 
(Il99ll) also concluded that the flu x density variatio ns can best 
be modelled with very hard spectra. [Gear et al.l (1 19941) also found 
that submillimetre spectra in a sample of 48 sources tended to be 
flat. Indeed, the distribution they found (their figure 7) is quite 
similar to ours (Fig.0. 

As noted above, alternative explanations for the a^F distri- 
bution cannot be totally excluded. However, we consider sources 
such as 2251 + 158 (Figs. [T6l and ITTb to be compelling evidence 
for the existence of very hard non-standard electron spectra. The 
implications of this possibility are discussed below. 

If the original a±m really is around -0.2, the distribution 
shown in Fig. \5\ becomes easy to understand: in many slightly 
older HFI components we see the expected spectral steepening 
from -0.2 to -0.7 and beyond. One way to test this hypothesis is 
to try to estimate the ages of the highest- frequency spectral com- 
ponents, either from their turnover frequencies or from the total 
flux-density monitoring: the older the component, the lower the 
turnover frequency and the steeper the spectra. However, this 
analysis is beyond the scope of this paper and reliable conclu- 
sions require more Planck data than presently exist. 

An electron energy spectral index as hard as 1.5 has ma- 
jor implications for the acceleration mechanisms dominating in 
blazar shocks. Particle acceleration efficiency is strongly gov- 
erned by the compression ratio of the flow, which in turn is con- 
fined by well-known plasma physics, preventing the spectral in- 
dex from becoming arbitrarily small. For a relativistic shock the 
traditional first-order Fermi mechanism usually produces spec- 
tral indices of 2.2 or so (and not 2.5, which follows from the 
nonrelativistic theory), but not smaller. There are a few options, 
however, for bypassing the s « 2.2 limit. For example, power- 
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law spectra with s > 1 have been found in certain kinds of 
oblique shocks (e.g.. iKirk & Heavens! 1 19891) . as well as in par- 
allel shocks in the case of large -angle particle sca ttering due 
to very strong turbulence (e.g., Stec ker et al.l 120071) . or when 
the compression r atio felt by the particles is higher than tha t 
felt by the plasma (lEllison et al.lll99(AIVirtanen & Vainioll2005h . 
The last alternative is an especially interesting one, because 
the combination of low density and relatively strong magnetic 
field — conditions likely to be found in Poynting-flux-dominated 
AGN jets — quickly leads to non-negligible Alfven speeds and 
enhancement of the compression ratio felt by the particle s and, 
consequently, to very hard particle spectra (iTam mi 2008aJ3, and 
references therein). High Alfven speeds and scattering centres 
that are not frozen into the plasma also enable second-order 
Fermi acceleration, which can produce powe r-law spectra even 
with s < 1 in relativistic parallel shocks (iVirtanen & Vainiol 
2005) and c an work on time - scales comparable to the fastest 
blazar flares (ITammi & Duffvll2009l) . 

As a simple order-of-magnitude reality check we can esti- 
mate the conditions required by the standard first-order accel- 
eration in a parallel step shock in the presence of an increased 
scatterin g-centre compress ion ratio due to turbulenc e trans- 
mission (iVainio et~ai1l2003l) . Following the analysis of Tammi 
(2008b), we can estimate that in order to achieve s « 1.5 in a 
shock with Lorentz factor r = 10, magnetic field strength of the 
order of 0.02-2 Gauss (depending on the composition and the 
density of the plasma) is needed. This is in good agreement with 
the usual modelling parameters of these sources, although we 
emphasise that the analysis used here is a simple one and may 
be limited to parallel shocks and weak turbulence. 

8. Conclusions 

We have presented the averaged Planck ERCSC spectra together 
with supporting ground and satellite observations obtained si- 
multaneously with the Planck scans. The 104 SEDs, supple- 
mented with archival data, demonstrate the usefulness of Planck 
data in determining the SED peak frequencies and fluxes and in 
modelling the spectral energy distributions in greater detail. In 
particular, the data demonstrate that the synchrotron spectrum 
contains contributions from several physically distinct AGN 
components: the LF and HF spectra are rarely smooth, except at 
the highest radio frequencies where the source component with 
the highest turnover becomes optically thin. The LF spectral in- 
dices cluster around a = 0, also indicating the superposition of 
many components with different turnover frequencies. In mod- 
elling the synchrotron and inverse-Compton SEDs, a multicom- 
ponent approach is therefore necessary. While physical mod- 
elling is beyond the scope of this paper, which uses averaged, 
rather than single-epoch, Planck ERCSC flux densities, we have 
shown some examples of multicomponent decomposition of the 
SEDs and the conclusions that can be drawn from the data. 

The sources display remarkable variability, which must also 
be taken into account in modelling the SEDs. We have presented 
some examples of total flux density monitoring to demonstrate 
how sources have been observed during Planck scans in different 
stages of activity. 

The Planck HFI spectra are remarkably flat, with only a mi- 
nority of sources having spectral indices steeper than the "canon- 
ical" a = -0.7, and many apparently having a t hi n around -0.2 
to -0.4. We suggest that the most likely interpretation for the 
very flat and straight synchrotron spectra observed by Planck is 
a very hard original electron spectral energy index, clearly be- 
low s = 2. While this possiblity must be confirmed by future 



data, an index s = 1.5 is compatible with the HF spectral index 
distribution seen in our sample. In a number of sources, energy 
losses steepen the HF spectra by Aa = 0.5 from about -0.3 to 
about -0.8, and further to -1 and beyond, which can explain 
the observed distribution. The hardness of the electron spectrum 
would also mean that the synchrotron SED peak is in most cases 
related to the maximum electron Lorentz factor achieved in the 
acceleration process. The hard spectra also require a rethinking 
of the acceleration processes dominating in the relativistic jets 
and shocks. 

For some sources the HF spectra are flat or even rising, in- 
dicating the presence of multiple synchrotron-emitting compo- 
nents peaking in the gigahertz to terahertz regime, and ongo- 
ing flaring at still higher frequencies. However, in many cases 
the Planck HFI spectrum appears to join the optical spec- 
trum smoothly, leaving little room for additional still higher- 
frequency components. In these cases, the inverse-Compton 
gamma-rays must originate in the synchrotron (i.e., shock) com- 
ponents seen in the Planck data. 
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Table 1. The complete 1 Jy northern sample of AGN. Column 8 refers to the 37 GHz observations of Metsahovi Radio Observatory. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Name 


J2000 Name 


Alias 


R.A.(J2000) 


Dec(J2000) 


1st scan 


2nd scan 


Save Uy] 


log 


Syn quality 


log vie 


IC quality 


0003-066 


0006-0623 


NRA05 


00:06:13.90 


-06:23:36.00 


2009-12-12 




2.06 


13.1 


A 






0007+106 


0010+1058 


IIIZW2 


00:10:31.01 


10:58:29.50 


2009-12-25 




1.04 










0048-097 


0050-0929 


PKS0048-097 


00:50:41.20 


-09:29:06.00 


2009-12-23 




1.34 


14.1 


A 






0059+581 


0102+5824 


TXS 0059+581 


01:02:45.76 


58:24:11.14 


2010-01-28 




2.97 


13.6 


B 






0106+013 


0108+0135 


OC012 


01:08:38.77 


01:35:00.30 


2010-01-04 




1.85 


12.7 


B 






J0125-0005 


0125-0005 




01:25:28.84 


-00:05:56.00 


2010-01-08 




1.24 


13.6 


B 






0133+476 


0136+4751 


DA55 


01:36:58.59 


47:51:29.10 


2010-01-29 




3.02 


12.8 


A 






0149+218 


0152+2207 




01:52:18.06 


22:07:07.70 


2010-01-24 




1.06 










0202+149 


0204+1514 


4C15.05 


02:04:50.41 


15:14:11.00 


2010-01-24 




1.64 










0212+735 


0217+7349 




02:17:30.81 


73:49:32.60 


2009-09-13 


2010-02-11 


2.05 


12.4 


B 






0224+671 


0228+6721 




02:28:50.05 


67:21:03.00 


2009-09-13 


2010-02-10 


1.28 










0234+285 


0237+2848 


4C28.07 


02:37:52.41 


28:48:09.00 


2010-02-03 




2.86 


13.1 


A 


22.1 


B 


0235+164 


0238+1636 




02:38:38.80 


16:36:59.00 


2010-02-01 




2.33 


13.1 


A 






0238-084 


0241-0815 




02:41:04.80 


-08:15:20.75 


2010-01-24 




1.33 










0306+102 


0309+1029 


PKS0306+102 


03:09:03.60 


10:29:16.00 


2010-02-06 




1.01 


12.9 


B 






0316+413 


0319+4130 


3C84 


03:19:48.16 


41:30:42.10 


2010-02-13 




14.01 






22.9 


B 


0333+321 


0336+3218 


NRAO140 


03:36:30.11 


32:18:29.30 


2010-02-15 




1.49 










0336-019 


0339-0146 


CTA026 


03:39:30.94 


-01:46:36.00 


2010-02-10 




2.21 


13.1 


B 






0355+508 


0359+5057 


NRAO150 


03:59:29.75 


50:57:50.20 


2009-09-13 


2010-02-19 


5.92 










0415+379 


0418+3801 


3C111 


04:18:21.28 


38:01:35.80 


2010-02-22 




5.97 










0420-014 


0423-0120 


OA129 


04:23:15.80 


-01:20:33.10 


2010-02-21 




7.27 


13.2 


B 






0430+052 


0433+0521 


3C120 


04:33:11.10 


05:21:15.60 


2010-02-24 




2.65 










0446+112 


0449+1121 


PKS0446+112 


04:49:07.67 


11:21:28.00 


2010-02-27 




1.41 










0458-020 


0501-0159 


PKS0458-020 


05:01:12.81 


-01:59:14.00 


2010-03-01 




1.86 










0507+179 


0510+1800 




05:10:02.37 


18:00:41.58 


2009-09-13 


2010-03-03 


1.14 










0528+134 


0530+1331 


PKS0528+134 


05:30:56.42 


13:31:55.15 


2009-09-14 


2010-03-07 


4.55 






21.8 


B 


0552+398 


0555+3948 


DA193 


05:55:30.81 


39:48:49.17 


2009-09-22 


2010-03-09 


3.18 






21.4 


B 


0605-085 


0607-0834 


PKS0605-085 


06:07:59.70 


-08:34:50.00 


2009-09-19 


2010-03-17 


1.73 


13.2 


B 






0642+449 


0646+4451 


OH471 


06:46:32.03 


44:51:16.59 


2009-09-30 


2010-03-16 


2.36 


12.5 


B 


22.7 


B 


0716+714 


0721+7120 




07:21:53.30 


71:20:36.00 


2009-10-04 


2010-03-10 


2.12 


14.3 


A 


23.4 


B 


0723-008 


0725-0054 


PKS0723-008 


07:25:50.70 


-00:54:56.00 


2009-10-09 


2010-04-02 


1.05 


13.0 


A 






0735+178 


0738+1742 


PKS0735+17 


07:38:07.40 


17:42:19.00 


2009-10-10 


2010-03-31 


1.93 










0736+017 


0739+0137 




07:39:18.03 


01:37:04.60 


2009-10-12 


2010-04-04 


1.85 


14.0 


B 






0748+126 


0750+1231 




07:50:52.05 


12:31:04.83 


2009-10-14 


2010-04-04 


2.42 


12.7 


A 


21.5 


A 


0754+100 


0757+0956 


OI090.4 


07:57:06.64 


09:56:34.90 


2009-10-15 


2010-04-06 


1.33 


13.0 


B 






0805-077 


0808-0751 




08:08:15.54 


-07:51:09.89 


2009-10-22 


2010-04-12 


1.05 


12.7 


B 






0804+499 


0808+4950 




08:08:39.67 


49:50:36.50 


2009-10-12 


2010-03-27 


1.30 


13.3 


B 






0823+033 


0825+0309 


PKS0823+033 


08:25:50.30 


03:09:24.00 


2009-10-24 


2010-04-14 


1.48 


13.7 


B 






0827+243 


0830+2411 


OJ248 


08:30:52.08 


24:10:60.00 


2009-10-21 


2010-04-09 


1.48 










0836+710 


0841+7053 


4C71.07 


08:41:24.37 


70:53:42.20 


2009-10-12 


2010-03-17 


1.80 


12.9 


B 


20.5 


B 


0851+202 


0854+2006 


OJ287 


08:54:48.80 


20:06:30.00 


2009-10-27 


2010-04-15 


3.90 


13.1 


B 


21.8 


B 


0906+430 


0909+4253 


3C216 


09:09:33.50 


42:53:46.08 


2009-10-23 


2010-04-09 


1.02 


13.2 


B 






0917+449 


0920+4441 




09:20:58.46 


44:41:53.99 


2009-10-25 


2010-04-11 


1.27 


12.7 


A 


21.7 


A 


0923+392 


0927+3902 


4C39.25 


09:27:03.01 


39:02:20.90 


2009-10-28 


2010-04-15 


7.06 






22.7 


B 


0945+408 


0948+4039 


4C40.24 


09:48:55.34 


40:39:44.60 


2009-10-31 


2010-04-18 


1.32 


13.1 


B 






0953+254 


0956+2515 




09:56:49.88 


25:15:16.05 


2009-11-08 


2010-04-27 


1.13 


13.1 


B 






0954+658 


0958+6533 


S40954+65 


09:58:47.20 


65:33:54.00 


2009-10-21 


2010-03-30 


1.03 


12.8 


B 


21.7 


B 


1036+054 


1038+0512 




10:38:46.78 


05:12:29.09 


2009-11-29 


2010-05-16 


1.13 











Table 1. continued. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ID 


(12) 


Name 


J2000 Name 


Alias 


R.A.(J2000) 


Dec(J2000) 


1st scan 


2nd scan 


Save [Jy] 




Syn quality 


log Vic 


IC quality 


TEX1040+244 


1043+2408 




10:43:09.00 


24:08:35.00 


2009-11-19 


2010-05-10 


1.02 










1055+018 


1058+0133 


OL093 


10:58:29.61 


01:33:58.80 


2009-12-07 


2010-05-24 


3.96 


12.7 


B 


22.5 


A 


J1130+3815 


1130+3815 


QSOB1128+385? 


11:30:53.28 


38:15:18.50 


2009-11-21 


2010-05-15 


1.16 


12.9 


B 






1150+812 


1153+8058 




11:53:12.50 


80:58:29.15 


2009-10-16 


2010-03-07 


1.11 










1150+497 


1153+4931 




11:53:24.47 


49:31:08.83 


2009-11-16 


2010-05-11 


1.40 










1156+295 


1159+2914 


4C29.45 


11:59:31.83 


29:14:44.00 


2009-12-04 


2010-05-31 


2.00 


13.4 


A 


23.5 


B 


1219+044 


1222+0413 




12:22:22.55 


04:13:15.78 


2009-12-29 




1.39 










1222+216 


1224+2122 


PKS1222+216 


12:24:54.51 


21:22:47.00 


2009-12-18 




1.16 


14.2 


B 


22.7 


B 


1226+023 


1229+0203 


3C273 


12:29:06.69 


02:03:08.60 


2010-01-02 




24.37 


14.1 


B 


20.8 


A 


1228+126 


1230+1223 


3C274 


12:30:49.42 


12:23:28.00 


2009-12-26 




14.04 


13.4 


B 






1253-055 


1256-0547 


3C279 


12:56:11.17 


-05:47:21.50 


2010-01-13 




17.20 


12.8 


B 


22.7 


B 


1308+326 


1310+3220 


AUCVn 


13:10:28.66 


32:20:43.80 


2009-12-21 




2.26 


12.6 


B 






1324+224 


1327+2210 




13:27:00.86 


22:10:50.16 


2010-01-05 




1.15 


12.8 


B 






1413+135 


1415+1320 




14:15:58.80 


13:20:24.00 


2010-01-25 




1.51 










1418+546 


1419+5423 


OQ530 


14:19:46.60 


54:23:14.00 


2009-12-06 




1.16 






20.6 


B 


1502+106 


1504+1029 


OR103 


15:04:24.98 


10:29:39.00 


2010-02-08 




1.30 










1510-089 


1512-0905 


PKS1510-089 


15:12:50.53 


-09:05:59.00 


2010-02-14 




2.48 


13.4 


B 






1546+027 


1549+0237 




15:49:29.44 


02:37:01.20 


2010-02-21 




2.42 


13.8 


B 






1548+056 


1550+0527 




15:50:35.27 


05:27:10.45 


2010-02-21 




2.04 










1606+106 


1608+1029 


4C10.45 


16:08:46.20 


10:29:07.80 


2010-02-24 




1.65 










1611+343 


1613+3412 


DA406 


16:13:41.00 


34:12:48.00 


2010-02-23 




3.03 


13.5 


B 






1633+382 


1635+3808 


4C38.41 


16:35:15.49 


38:08:04.50 


2010-03-02 




2.90 






21.6 


B 


1638+398 


1640+3946 




16:40:29.63 


39:46:46.03 


2010-03-04 




1.18 










1642+690 


1642+6856 




16:42:07.85 


68:56:39.76 


2009-11-06 


2010-01-17 


1.74 










1641+399 


1642+3948 


3C345 


16:42:58.81 


39:48:37.00 


2010-03-05 




7.89 


12.2 


A 






1652+398 


1653+3945 


MARK501 


16:53:52.20 


39:45:36.00 


2010-03-09 




1.00 










1739+522 


1740+5211 


S41739+52 


17:40:36.98 


52:11:43.00 


2010-03-30 




1.45 










1741-038 


1743-0350 


PKS1741-038 


17:43:58.86 


-03:50:04.60 


2009-09-13 


2010-03-18 


4.39 










1749+096 


1751+0939 


PKS 1749+096 


17:51:32.70 


09:39:01.00 


2009-09-13 


2010-03-21 


4.50 


13.0 


A 


23.2 


A 


1803+784 


1800+7828 


S51803+784 


18:00:45.40 


78:28:04.00 


2009-10-11 


2010-02-02 


1.90 


13.5 


A 


23.1 


A 


1807+698 


1806+6949 


3C371.0 


18:06:50.70 


69:49:28.00 


2009-10-21 


2010-01-07 


1.37 


14.3 


A 


23.1 


B 


1823+568 


1824+5651 


4C56.27 


18:24:07.07 


56:51:01.50 


2010-04-17 




1.60 


13.0 


B 






1828+487 


1829+4844 


3C380 


18:29:31.80 


48:44:46.62 


2010-04-11 




2.38 


14.0 


B 






J184915+67064 


1849+6705 




18:49:15.89 


67:06:40.90 


2009-10-12 


2009-12-25 


1.82 


12.6 


A 






1928+738 


1927+7358 


4C73.18 


19:27:48.50 


73:58:01.60 


2009-09-28 


2010-01-18 


2.61 






21.0 


B 


1954+513 


1955+5131 




19:55:42.74 


51:31:48.55 


2009-10-09 


2010-05-06 


1.35 


13.7 


B 






2007+776 


2005+7752 


S52007+77 


20:05:31.10 


77:52:43.00 


2009-09-25 


2010-01-27 


1.40 


13.1 


B 






2005+403 


2007+4029 




20:07:44.94 


40:29:48.60 


2009-10-17 


2010-05-02 


2.04 










2021+614 


2022+6136 


OW637 


20:22:06.68 


61:36:58.80 


2009-12-20 


2010-05-26 


1.25 










2037+511 


2038+5119 




20:38:37.04 


51:19:12.66 


2009-11-25 


2010-05-18 


2.46 


13.2 


B 






2121+053 


2123+0535 




21:23:44.52 


05:35:22.09 


2009-11-02 


2010-05-08 


1.74 










2131-021 


2134-0153 


PKS213 1-021 


21:34:10.31 


-01:53:17.24 


2009-11-02 


2010-05-09 


1.46 










2134+004 


2136+0041 


OX057 


21:36:38.59 


00:41:54.21 


2009-11-04 


2010-05-10 


3.34 






21.0 


A 


2136+141 


2139+1423 




21:39:01.31 


14:23:36.00 


2009-11-11 


2010-05-16 


1.29 


12.5 


B 






2145+067 


2148+0657 




21:48:05.46 


06:57:38.60 


2009-11-10 


2010-05-16 


7.37 










2200+420 


2202+4216 


BLLAC 


22:02:43.30 


42:16:39.00 


2009-12-13 


2010-06-06 


3.35 


13.8 


B 


22.6 


A 


2201+315 


2203+3145 


4C31.63 


22:03:14.98 


31:45:38.30 


2009-12-02 


2010-05-30 


2.53 


14.2 


B 


20.7 


B 


2201+171 


2203+1725 




22:03:26.89 


17:25:48.20 


2009-11-21 


2010-05-24 


1.21 


13.5 


B 
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(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ID 


(12) 


Name 


J2000 Name 


Alias 


R.A.(J2000) 


Dec(J2000) 


1st scan 


2nd scan 


Save [Jy] 


logv, 


Syn quality 


log Vic 


IC quality 


2216-038 


2218-0335 




22:18:52.04 


-03:35:36.90 


2009-11-14 


2010-05-21 


1.41 


14.0 


B 






2223-052 


2225-0457 


3C446 


22:25:47.26 


-04:57:01.40 


2009-11-15 


2010-05-23 


5.58 


12.7 


B 


22.0 


A 


2227-088 


2229-0832 




22:29:40.09 


-08:32:54.50 


2009-11-14 


2010-05-23 


1.94 






20.7 


B 


2230+114 


2232+1143 


CTA102 


22:32:36.41 


11:43:50.90 


2009-11-27 


2010-05-31 


3.84 


13.2 


B 


22.6 


B 


2234+282 


2236+2828 




22:36:22.47 


28:28:57.41 


2009-12-10 




1.03 


13.0 


B 






2251+158 


2253+1608 


3C454.3 


22:53:57.75 


16:08:53.60 


2009-12-06 




9.42 


13.9 


B 


22.6 


A 


4C45.51 


2354+4553 




23:54:21.68 


45:53:04.00 


2010-01-11 




1.12 










2353+816 


2356+8152 




23:56:22.79 


81:52:52.26 


2009-09-19 


2010-02-08 


1.59 
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Table 4. Simultaneity of high energy data to Planck observations displayed in the SEDs. 
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z /m 


UZZ4+0 / 1 
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sim_lst 
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• 

sim_lst 


• 

sim_lst 


sim_lst 
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sim_lst 


sim_lst 


sim_lst 


nooo HQ/1 
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z /m 


U333+3Z1 


U33o+3zlo 
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U33y— U140 






sim_lst 


U3jj+jUo 


U3 0^+ JUJ / 
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Table 4. continued. 
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Fig. 18. Left-panel: The SED of the source 0003-066. Grey circles show the historical data. The red circles show data simultaneous 
to the Planck observations. The dotted and dashed lines show the second and third degree polynomials, respectively, fitted to the 
synchrotron and IC bumps in the SED. Right-panel: The radio spectrum of 0003-066. Red circles, LF data simultaneous to Planet, 
red stars, ERCSC LFI data; blue circles, HF data simultaneous to Planck; blue stars, ERCSC HFI data. The dashed and dotted lines 
are fits to simultaneous LF and HF data. 
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Fig. 32. 0306+102 
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Fig. 38. 0420-014 
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Fig. 42. 0507+179 
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Fig. 43. 0528+134 
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Fig. 44. 0552+398 
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Fig. 46. 0642+449 
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Fig. 47. 0716+714 
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Fig. 50. 0736+017 
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Fig. 53. 0805-077 
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Fig. 56. 0827+243 
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